Can one identify the intrinsic structure of the yrast states in 48 Cr after the 

backbending? 
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The backbending phenomenon in 48 Cr has been investigated using the recently developed Pro- 
jected Configuration Interaction (PCI) method, in which the deformed intrinsic states are directly 
associated with shell model (SM) wavefunctions. Two previous explanations, (i) K = band cross- 
ing, and (ii) K = 2 band crossing have been reinvestigated using PCI, and it was found that both 
explanations can successfully reproduce the experimental backbending. The PCI wavefunctions in 
the pictures of K = band crossing and K = 2 band crossing are highly overlapped. We conclude 
that there are no unique intrinsic states associated with the yrast states after backbending in 48 Cr. 

PACS numbers: 21.60.Ev, 21.60.Cs, 21.10.Re, 27.40,+z 



I. INTRODUCTION 

The backbending of 48 Cr has been observed more than 
10 years ago [J, HI, but its interpretation remains contro- 
versial and challenges the existing nuclear models. Shell 
Model (SM) calculations have reproduced very well the 
yrast states of 48 Cr [1,0], but it is difficult to provide the 
physical insight because the laboratory frame wavcfunc- 
tion doesn't imply any information associated with the 
deformed intrinsic structure. The cranked Hartree Fock 
Bogoliubov(CHFB) method is a complementary theory 
Q often used to analyze the deformed intrinsic states. 
According to one CHFB analysis [4], 48 Cr is an axial rotor 
up to the backbending, after that the system changes to a 
spherical shape. An alternative and more detailed CHFB 
analysis @ shows that the backbending of the 48 Cr is not 
associated with the single particle level crossing and that 
the intrinsic configuration remains unchanged. 

The Projected Shell Model (PSM) @, H is an alterna- 
tive technique that mixes the best intrinsic shell model 
configuration with other associated particle-hole configu- 
rations. A PSM analysis 9] indicates that the backbend- 
ing in 48 Cr is due to a band crossing involving an excited 
4-quasiparticle (qp) band with K = 2, which represents 
a configuration of broken neutron and proton pairs. Un- 
fortunately, such an explanation is not confirmed by re- 
cent SM [1 03 and CHFB 0, Q analyses. In addition, 
the same PSM Ref. uses the Generator Coordinate 
Method (GCM) that provides the picture of a spheri- 
cal band crossing, which in some sense is similar to the 
result of CHFB. Obviously, the above apparently con- 
flicting explanations coming from different models need 
to be reinvestigated. 

In the present paper, the backbending in Cr is inves- 
tigated in the framework of the newly developed method 
called Projected Configuration Interaction (PCI) fill. [l2l|. 
The PCI basis is built from a set of Slater determinants 
(SD). Those SDs may have different shapes, including the 
spherical shape. Hence, the nuclear states with different 
intrinsic shapes can be mixed by the residual interaction. 



By using the same SM Hamiltonian, the PCI results can 
be directly compared with those of full Shell Model cal- 
culations. Moreover, PCI uses deformed single particle 
bases and therefore, the physics insight of the results can 
be clearly analyzed. Different PCI bases can be built 
in such a way that they reflect the nature of the intrin- 
sics states found in previous studies, such as CHFB (or 
GCM) and PSM. PCI wavefunctions were shown to be 
very good approximations to those of full SM, and they 
can be obtained using different bases. Thereafter, over- 
laps among PCI wavefunctions can be calculated and an- 
alyzed to determine the validity of various explanations. 
These features suggest that PCI could shed new light 
on the interesting phenomenon of backbending in 48 Cr. 
Other models using similar techniques includes the family 
of VAMPIR [ll|, and the quantum Monte Carlo diago- 
nalization (QMCD) method [T3|. 

The paper is organized as follows: Section II provides 
a short introduction into the basics concepts used by the 
PCI method that would be required for an understand- 
ing of the arguments used later in the analysis. Section 
III is devoted to the analysis of the contribution of dif- 
ferent intrinsics set of states to the backbending in 48 Cr. 
Section IV summarizes the conclusion of our study. 



II. BASIC CONCEPTS OF THE PCI METHOD 

For completeness we give here a brief introduction of 
the PCI method (see Refs. [ill, EH f° r more details). 
The deformed single particle (s.p.) states need to be 
generated either by HF calculation or from a deformed 
s.p. Hamiltonian. For simplicity, here we take the latter 
approach, and the s.p. Hamiltonian can be written as 



(1) 



where h sp ^ = eic\ci is the spherical s.p. Hamilto- 
nian assumed to have the same eigenfunctions as the 
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spherical harmonic oscillator, energies are properly ad- 
justed such that the SD with lowest energy is close to 
the HF vacuum. For the p/-shell we use e/ 7 , 2 = O.OMeV, 
e P3/2 = 4.5MeV, e fs/2 = 5.0MeV and e Pl/2 = 6.0MeV. 
In Eq. (HJ) £2, £4 are the quadrupole and hexadecapole 
deformation parameters, Pi are Legendre polynomials, 
p = r/b, and we take b = 1.01A 1 / 6 for the harmonic 
oscillator parameter [3, 0] . 

The Slater determinants can be built with deformed 
s.p. states. Following our previous papers [HI, the 
general structure of the PCI basis can be written as 
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where |ftj,0) (i = 1, ...N) is an optimal [12] set of start- 
ing states having different deformations. Assuming that 
these \k, 0) are found (in what follows we skip the sub- 
script i to keep notation short), a number of relative np- 
nh SDs, \k, j), on top of each \k, 0) are added to the SD 
basis selected the constraint PJJ 



AE 
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Ej) 2 +4\V\ 2 ) >E cut . (3) 



Here E = (k,0\H\k,0), Ej = {K,j\H\n,j) and V = 
(k, 0\H\n,j). The PCI basis is then obtained by project- 
ing the selected SDs onto good angular momentum. The 
wavefunctions, as well as the energy levels, are obtained 
by solving the following generalized eigenvalue equation: 



{H I KHf -E I Nf M ,)fi f =Q. 



(4) 



Here Hl„, and N* , are given by 



= (k\hp£ k ,W), n. 



I 



= (*\Pkk>W), (5) 



where P— K , is the angular momentum projection op- 
erator, and H is the shell model Hamiltonian. In this 
study we take the KB3 interaction [15jj, which has been 
used by Caurier et al in their shell model calculations of 
48 Cr [4[. However, results similar to those reported here 
are provided by other interactions, such as KB3G [l6| or 
GXPF1A 171. 



III. THE PCI ANALYSIS OF THE 
BACKBENDING IN 48 CR 



To get some insight into the structure of the states 
contributing to the backbending in 48 Cr we chose two 
basic |k, 0) SDs. The first |/ci,0) is a K — configu- 
ration with all 8 valence nucleons occuping the |f2| = 
1/2(1/7/2) orbits and the |0| = 3/2(l/ 7 / 2 ) orbits, as 
shown in Fig. QJ,. The deformation of |«i,0) is given 
by 62 = 0.19 and 64 = —0.05. This deformation was 



TABLE I: PCI bases used for the backbending study in 48 Cr. 
K = and K — 2 configurations are shown in Fig. [T] 



Basis 



M) 



g.s. 



K €2 £4 
0.19 -0.05 
143 



A 



K €2 £4 

0.00 0.00 
192 



B 



K €2 £4 
2 0.00 0.00 
186 



C 



£2 

0.19 
149 



£4 

-0.05 




FIG. 1: (color online) (a) The K = and K = 2 configu- 
rations used in the present calculations of 48 Cr. All 4 levels 
are from I/7/2 subshell. (b) Potential energy curves of the 
K — configuration in (a) as functions of £2. The solid lines 
show the projected energies and the dashed line, unprojected 
energy. £4 was chosen to minimize the potential energy for 
each £2. KB3 interaction was used. 



obtained by determining the minimum of the energy sur- 
face of (K = 0|ffKB3|-^ = 0) as a function of e 2 and 64, 
(See the dashed line in Fig. [Tp). A minimum energy of 
(ki,0|.Hkb3|«i,0) = -28.296 MeV was found, which is 
close to the HF energy -28.423 MeV. |«i,0) is believed 
to be responsible for the low-spin yrast states before the 
backbending, and the corresponding PCI basis that in- 
cludes the np-nh states selected by Eq. ([3]) is denoted as 
'g.s.', and it is shown in Table HI 

The second basic SD, |k 2 ,0), is chosen to describe the 
high spin states in 48 Cr, after the backbending. Accord- 
ing to the previous studies mentioned in the introduction, 
there are at least two candidate configurations for \k2, 0). 
The first, suggested by the CHFB calculations, according 
to which the backbending in 48 Cr can be explained with- 
out a band crossing [(| , but by a shape changes from well 
deformation to spherical Q. Our PCI calculations seem 
to be in agreement with this interpretation. The pro- 
jected energy curves for the K = configuration at each 
even-spin are shown in Fig. [Tp (solid lines). One can 
see that the deformation at minimum decreases as the 
spin increases. Guided by this result, one can establish 
a possible |k2,0), which has the same configuration as 
«i,0), but whose shape is almost spherical. Such choice 
of «2,0) SD would be consistent with the GCM inter- 
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FIG. 2: (color online) (a) PCI energies E(I) as functions of 
spin for the bases listed in Table H] (b) Overlaps among the 
PCI wavefunctions for bases A, B, and C. 



pretation, in which HF vacua with different deformations 
are included in the basis. However, in PCI one can fur- 
ther include particle-hole excitations on top of each Ik, 0). 
On should recall that in the PSM interpretation [9[ the 
backbending in 48 Cr is caused by a K — 2 band crossing. 
Therefore, PCI seem to be ideally suited to include and 
analyze another possible |k2,0) with K = 2. Its config- 
uration is also shown in Fig. [TJi, and its structure can 
be chosen cither spherical (see column B in Table fl} or 
of the same deformation as that of 0) favored by the 
PSM approach (see column C in Table HJ. In an attempt 
to find the optimal structure of yrast states in 48 Cr after 
the backbending, we considered all three possibilities of 
|« 2 ,0), labeled with 'A', 'B' and 'C in TableU 

The np-nh SDs built on top of each \k,Q) are 

selected by setting E cut — 0.5 keV in Eq. @. Conse- 
quently, the number of selected j) is 142, and those 
of \n2,j) for A, B and C are 191, 185 and 148, respec- 
tively. Adding the \k, 0} itself, the total number of the 
selected SDs, N K , for each |k, 0} is listed in TableU] 

These PCI bases were used to calculate all even-spin 
states from / — to / = 16. The calculated lowest ener- 
gies with basis g.s., A, B and C as functions of spin are 
shown in Fig. Uta). The energies calculated with basis 
g.s. alone increases very smoothly with the spin, with- 
out exhibiting any backbending. For the bases A, B and 
C, all the corresponding energy curves cross that of the 
g.s. band after spin / = 10. To get a better understand- 
ing of the relation among those bases, we calculated the 
overlaps between corresponding wavefunctions. Let's de- 



FIG. 3: (color online) (a) E2 transition energies E(I) — E(I — 
2) vs spin I. Experimental data are taken from Ref. [2J, and 
the result of the Full SM from Ref. 4]. (b) Overlaps among 
the PCI wavefunctions with basis sets g.s.+A, g.s.+B and 
g.s.+C. 

note the PCI wavefunctions as \ty K (I)), where n refers 
to certain basis or combination. For instance, k can be 
'g.s.', 'A', or 'g.s.+A', etc. Overlaps of \{V A (I)\V B (I))\, 
\(y A (I)\$c(I))\ and \{Vb(I)\Vc(I))\ as functions of/ 
are plotted in Fig. &b). The surprising result is that 
these overlaps are unexpectedly very large, which means 
quite different intrinsic bases would generate almost the 
same wavefunctions after the angular momentum projec- 
tion. In particular, one should note the large overlap 
|(*a(-0I*b(-0)| « 1 at 7 > 10, while without the an- 
gular momentum projection bases A and B are strictly 
orthogonal due to different K values . 

The E 1 (I) = E(I) - E(I - 2) energies calculated us- 
ing combinations of bases, g.s.+A, g.s.+B and g.s.+C 
are shown in Fig. [3^a). The backbending phenomenon 
seems to be easily reproduced by all those bases. As al- 
ready discussed, the g.s.+A basis is qualitatively similar 
to that used in the CHFB and GCM analyses, while the 
g.s.+C basis follows the scenario proposed by the PSM 
analysis. Furthermore, as shown in Fig. OUb), the over- 
laps \{^ g. s .+A(i)\^ g.s.+c(I))\ are very large, at least 92%. 
Therefore, we come to the conclusion that the appar- 
ently contradictory CHFB and PSM explanations, actu- 
ally look qualitatively equivalent one to each other. One 
should note, however, that the backbending described 
with the g.s.+C basis is quantitatively not as good as 
that descibed with the g.s.+A basis. The reason could 
be the large deformation of C basis. Changing the defor- 
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FIG. 4: (color online) (a) Yrast state energies vs spin obtained 
with PCI (open symbols) and Full SM(solid symbols), (b) 
The BE2 values with the same wavefunctions as in (a). Result 
of the Full SM are taken from Ref. @|. 



mation of C basis to spherical, one gets the g.s.+B basis, 
and the result is improved. This feature supports the 
idea that the shape of 48 Cr reduces after backbending, 
as has been pointed out in Ref. [H, 0] . 

Notice that the results obtained with bases g.s.+A and 
g.s.+B arc almost identical, although the intrinsic bases 
A and B have a completely different structure. Our calcu- 
lations show that there are also other intrinsic bases with 
K different of and 2 that can reproduce the backbend- 
ing, and whose corresponding wavefunctions are almost 
equivalent to the basis g.s.+A or g.s.+B when projected 
on good angular momenta. In other words, one can not 
find a unique intrinsic state for the yrast states in Cr 
for I = 12 — 16. One can get some insight into the ap- 
parent irrelevance of the intrinsic structure at high spins 
by analyzing the case of / = 16, which is the band termi- 
nation state. In the space of Trlfy 2 vlfy 2 , there is only 
one SD that reaches the maximum K = 16, showing that 
only one / = 16 state can be constructed from that space. 
On the other hand many irlfy 2 



v\fy 2 SDs with various 
K values can be projected onto good /. Therefore, we 
have many projected states with / = 16, and our cal- 
culations prove that they are exactly identical and the 
projected energy is —18.342 MeV with KB3 interaction. 



Therefore, one can use any one of the wlfy 2 



1/1/4/2 SDs 



to reproduce the state of band termination. 

Fig. [3] shows good comparison of the E2 transition 
energies obtained with the particular choice of bases de- 
scribed above and those of the full SM calculations. One 
would like to further compare and validate the structure 
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FIG. 5: (color online) Overlaps O(I) in Eq[6] between the 
projected intrinsic configurations with K = and K = 2 
shown in Fig[l}i for spherical and deformed (e? = 0.19 and 
£4 = —0.05) cases. 



of different approaches with the full SM results. It is dif- 
ficult to calculate the overlap between PCI and full SM 
wavefunctions because the large number of norm matrix 
elements, Apci x Asm that has to be calculated (see Eq. 
([5])). Here Npci and Asm are dimensions of the PCI and 
full (M-scheme) SM spaces, where Asm is a very large 
number. However, the PCI approximation can be fur- 
ther checked by comparing the energies as well as the 
BE2 transition probabilities with those of the full SM 
calculations. The results are shown in Fig. |4] The PCI 
energies lie 500 keV on average above the exact values, 
but they were calculated using less than 400 SDs (see Ta- 
ble. HI, while in full SM calculations the number of the 
SDs used is 1963461. For the B(E2) calculations we used 
the 'canonical' effective charges, i.e. 1.5e for protons and 
0.5e for neutrons, which are the same as those used in 
Refs. [3, 4]. The BE2 values obtained with PCI are very 
close to the exact ones, except those for I = 6, 8 and 
10, which are slightly larger. That is very likely due to 
the fixed deformation we used for the selected 'g.s' basis. 
As shown in Fig. [T}), the deformation reduces gradu- 
ally as the spin increases, conclusion also supported by 
the CHFB [| and the cranked Nilsson-Strutinsky (CNS) 
[HI calculations. Therefore we overestimated the BE2 
values at / = 6, 8 and 10. 

On a final note the intrinsic states play a key role in 
studying the physics of the nuclear system. Even if they 
are orthogonal, their overlap can be quite large after the 
angular momentum projection. Here is a simple exam- 
ple. The K = and K — 2 configurations in Fig. QJ, 
were projected on to good angular momentum assum- 
ing spherical and deformed (e 2 — 0.19 and 64 = —0.05) 
shapes, respectively. The overlap 



0(1) 



(K = 0\Pi 2 \K = 2) 



y/{K = 0\P^\K = 0){K = 2\Pi 2 \K = 2) 



(6) 



was calculated and plotted in Fig. [S] One can see 
that the deformed overlap O(I) is always smaller than 
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the spherical one, indicating that the deformed intrin- 
sic states can be more clearly defined than the spherical 
ones. In the low-spin region the overlaps O(I) are small, 
and the two intrinsic configurations can be easily distin- 
guished. But at high spin, 0{I) are much larger for both 
spherical and deformed cases, which means that the in- 
trinsic states can not be clearly identified, at least for the 
case of 48 Cr. 



backbending in 48 Cr has been reproduced by using var- 
ious PCI bases, carefully selected to reflect the physics 
of the intrinsic states found by the CHFB (GCM) and 
PSM analyses of this case. Using the PCI capabilities of 
mixing these bases we show for the first time that the 
backbending pictures proposed by the CHFB and PSM 
methods are qualitatively equivalent. Our analysis sup- 
ports the conclusion that there is no unique intrinsic state 
for spins larger than 10 in 48 Cr. 



IV. SUMMARY 



In summary, the backbending in 8 Cr has been stud- 
ied with a recently developed Projected Configuration 
Interaction (PCI) method. PCI uses the same realistic 
Hamiltonians and valence spaces as the SM calculations, 
but only a set of properly selected SDs with different 
deformations and associated np-nh configurations. The 
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